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Abstract 

We present an analytic three-loop calculation for thermodynamic quantities of the 
0(n) symmetric 0^ theory below Tc within the minimal subtraction scheme at fixed 
dimension d = 3. Goldstone singularities arising at an intermediate stage in the 
calculation of 0(n) symmetric quantities cancel among themselves leaving a finite 
result in the limit of zero external field. Prom the free energy we calculate the three- 
loop terms of the amplitude functions fip, Fj^ and F_ of the order parameter and the 
specific heat above and below Tc, respectively, without using the e = 4— d expansion. 
A Borel resummation for the case n = 2 yields resummed amplitude functions /(^ and 
F_ that are slightly larger than the one-loop results. Accurate knowledge of these 
functions is needed for testing the renormalization-group prediction of critical-point 
universality along the A-line of superfluid ^He. Combining the three-loop result for 
F- with a recent five-loop calculation of the additive renormalization constant of 
the specific heat yields excellent agreement between the calculated and measured 
universal amplitude ratio /A~ of the specific heat of ^He. In addition we use 
our result for to calculate the universal combination Rc of the amplitudes of 
the order parameter, the susceptibility and the specific heat for n = 2 and n = 3. 
Our Borel-resummed three-loop result for Rc is significantly more accurate than 
the previous result obtained from the e-expansion up to O(e^). 
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1 Introduction 



Field-theoretic perturbative calculations of the critical behavior of 0(n) sym- 
metric systems with n > 1 below are known to be considerably more 
complicated than those of Ising-like (n = 1) systems. The difficulties for n > 1 
are due to the existence of transverse fluctuations of the order parameter (in 
addition to the ordinary longitudinal fluctuations of n = 1 systems) which, 
in the long-wavelength limit and at vanishing external field, have a vanishing 
restoring force. This implies the existence of massless Goldstone modes [1,2] 
which yield infinite transverse and longitudinal susceptibilities [3] and which 
cause infrared (Goldstone) singularities at intermediate stages of perturbative 
calculations of all thermodynamic quantities on the coexistence curve below 
Tc- These complications for n > 1 have prevented calculations of the equation 
of state and of amphtude functions below Tc to higher than two-loop order [4- 
6] , in contrast to the case n = 1 below Tc where accurate Borel resummation 
results based on five-loop perturbation theory are available [7-10]. 

Higher-order calculations of the amplitude functions below Tc for the case 
n — 2 are of primary importance in view of the proposed theoretical research 
[11] parallel to the considerable experimental effort [12] to test the universality 
prediction of the renormalization-group theory [13] along the A-line of *^He. 
The amplitude functions contain the information about universal ratios [13] 
of leading and subleading amplitudes near criticality that have previously 
been used as fit parameters in the data analysis [14,15] because of the lack 
of accurate predictions for these universal ratios. The amplitude functions are 
also a crucial ingredient in a nonlinear renormalization-group analysis [16-18] 
(equivalent to a resummation of the whole Wegner correction [19] series) in a 
wide temperature range including non-asymptotic and non-universal effects. 

In the present paper we perform the next necessary step towards the goal of 
an accurate determination of amplitude functions by presenting the analytic 
results of a three-loop calculation within the 0(n) symmetric 0*^ theory above 
and below for general n. Specifically, from the free energy we shall derive the 
amplitude function of the order parameter and the amplitude functions F± 
of the specific heat above and below in the limit of vanishing external field. 
The conceptual framework of our calculation is the minimally renormalized 
massive 0*^ field theory at fixed dimension d = 3 [5,17,20], without using the 
e — 4 — d expansion, involving an appropriately defined pseudo-correlation 
legth ^_ [5,6] that is applicable to both n — 1 and n > 1 below Tp. Our 
d — 3 approach differs from Parisi's [21] original suggestion of the d = 3 field 
theory and from the d = 3 field theory of subsequent work [7,9,22-24] where 
renormalization conditions rather than the minimal renormalization scheme 
were used. Our combination of the minimal subtraction scheme with the field 
theory at fixed 2 < d < 4 has several advantages, such as the independence 
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of renormalization constants on whether T > Tc or T < Tc which imphes a 
natural decomposition of correlation functions into amplitude functions and 
exponential parts where the structure of the latter is independent of whether 
T > Tc or T < Tc. The same exponential parts can then be used, without 
modification, in extensions of the theory to critical dynamics [25] or to critical 
phenomena in confined systems above and below Tc [26]. 

Our approach has been successfully employed recently [6] in deriving various 
amplitude functions, including and F±, for general n below Tc in two- loop 
order. In three-loop order, the technical difficulties related to the removal of 
ultraviolet divergences in three dimensions and to the treatment of spurious 
infrared (Goldstone) divergences are substantially greater. The concept of the 
minimally renormalized d = 3 field theory turns out to constitute an appro- 
priate framework for coping with these difficulties. This requires to further 
develop new integration techniques based on recent advances by Rajantie [27] 
in the analytic evaluation of thrcc-loop integrals in three dimensions. In par- 
ticular we succeed in calculating two "Mercedes" diagrams with two different 
(longitudinal and transverse) masses that were not considered previously. 

The perturbation series of the 0^ theory are known to be divergent and to 
require resummations [28] in order to yield reliable results. The three-loop 
terms derived in the present paper make possible to perform Borel resumma- 
tions that lead to results with reasonably small error bars. In a related paper 
[10] this was shown for F_. In the present paper we demonstrate this for the 
case of f^. For the example n — 2 the Borel resummed three-loop amplitude 
functions F_ and of the specific heat below Tc and of the order param- 
eter turn out to be slightly larger than the one-loop results. Application of 
the Borel resummed amplitude function F_ to the universal amplitude ratio 
A^/A~ of the asymptotic specific heat yields the theoretical prediction [10] 
A^/A~ — 1.056±0.004 which is in excellent agreement with the high-precision 
experimental result [15] A+fA' = 1.054 ± 0.001 for ^He near the superfiuid 
transition obtained from a recent experiment in space. Furthermore we use 
our result for to calculate the universal combination Rc [13] of the ampli- 
tudes of the order parameter, the susceptibility and the specific heat for n = 2 
and n — 3. Our Borel resummed three- loop result for Rc is significantly more 
accurate than the previous result obtained from the e-expansion up to O(e^) 
[29]. 



2 Bcire Helmholtz free energy 

Throughout this paper we shall use the notation and definitions of Ref . [6] . We 
start from the 0{n) symmetric 0^ model with the Landau-Ginzburg- Wilson 



3 



functional 

n{U^)} = j d'^x Qro0^ + i E(V0o^)' + uMf - ho ■ 0o) (1) 

y \ i / 

for the n-component field 0o(x) = (0oi(x), . . . , 0on(x)) in the presence of the 
homogeneous external field ho = {ho, 0, . . . , 0). The spatial variations of 0o(x) 
are restricted to wavenumbers less than some cutoff A. We are interested in 
the bulk Helmholtz free energy ro{ro,Uo, Mq, A) per unit volume with Mq = 
(0oi). To is obtained from the negative sum of all one-particle irreducible 
vacuum diagrams. The structure of the analytic expression is (apart from an 
unimportant additive constant) 

1 1 ? 

ro(ro, uo, Mo, A) = -roM^ + uoM^ + 2 J (""o^ + 

p 

A 

1 r 

+ -(n - 1) y In (for + p') + «o^i'^ (ro, uo, Mo, A) 
p 

+ ulxl,'\ro,Uo,Mo,A) + 0{ul) (2) 

where = (271)"^ ^ d'^ p means integration up to |p| < A. The terms 

'UoXQ^''(ro, mq, Mq, A) and ■UQXQ^''(ro, MQ) Mq, A) represent the two- and three- 
loop contributions shown in Fig. 1, with longitudinal and transverse propaga- 
tors Gl{p) = (^oL + P^)~^ and Gt{p) = (for + P^)'^ where 

foL = ro + l2uoMl , for = ro + AuoM^ (3) 

(compare Eq. (27) and Fig. 3 of Ref. [6] at /c = 0). For the integral expressions 
of Xo see Appendix A. 

We consider the limit A — > 00 (although neglecting cutoff effects may not be 
justified in certain cases [30-32]). The right-hand side of Eq. (2) is ultraviolet 
divergent for A — cxd at any c? > 2. To absorb these divergences requires 
two different steps: a shift of the temperature variable ro (mass shift) and 
a subtraction from the free energy. The first step is to turn to the shifted 
variable ro — roc where roc{uo,K) is the critical value of ro [20]. Substituting 
ro = ro — roc + Toduo, A) into Fo yields the function 

f o(ro - roc, Uo-, Mq, A) = ro(ro - roc + roc(Mo, A), mq, Mq, A) . (4) 

(Because of the non- analytic uq dependence of roc, the function f does not 
have an expansion in integer powers of Uo-) This function is still divergent for 
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UoX^ = .iuo 




XJ 



+ g(«-l)( 



+ 1728«oMo 




+ 3(«-l)t 



+ g('^-l) 



+ — (n - r 






+ g(«-l) 



+ —{n - 1) 



+ — (n - 1) 



Fig. 1. Two- and three-loop vacuum diagrams determining the Helmholtz free en- 
ergy ro(ro, Mo, Md, A), Eq. (2). The hues denote longitudinal and transverse prop- 
agators, Gl{jp) = = {foL+P^)~^ and Gt{p) = — \ — = (for The 

integral expressions are given in Appendix A. The analytic results for the last two 
"Mercedes" diagrams at c? = 3 are given in Eqs. (A. 7) and (A. 8). 

d > 2 in the limit A — > oo at fixed ro — roc, uq, Mq. We have verified up to 
three-loop order that these ultraviolet divergences have a regular dependence 

on Tq — roc (up to linear order in Tq — roc for d < 4 and up to quadratic 
order in Tq — Tqc for d = 4). Therefore it suffices to perform additional overall 
subtractions from Tq that are regular in Tq— roc- They do not affect the singular 
part of the temperature dependence of the free energy. The resulting singular 
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part of the bare free energy is finite in the hmit A — > oo for 2 < d < 4 at fixed 
^0 - ^Oc- 

Here we use the prescriptions of dimensional regularization at A = oo and 
employ the critical parameter tqc in the form [20,33] 

roc{uo,e) =uJ^S{e) (5) 



where the function S{e) is finite for e > except for poles at = 4 — 2//, / = 2, 
3, ... In our version of the d — 2) theory it is sufficient to work near d = 3 and, 
instead of Tq — roc, to use a shifted variable Tq as defined by Tq = Tq + 5ro(wo, 
with a simplified 5ro(Mo, e) that contains only the d — ?> pole of roc (and not 
the poles at ci; > 3 with / > 3). Thus we use [6-8] 



2/e 



5rQ{uQ, e) = — (n + 2) + C{n) Uq 



2A 



(6) 



with e — A — d and the finite constant 

n + 2 



An 



1 -CEuier + ln — -21n24 

y 



(7) 



(the final results for the amplitude function do not depend on the choice of 
C{n) [6,20]). Correspondingly, instead of Tql and for, we use the longitudinal 
and transverse parameters 

roL = + l'2uoMl , r^r = r^, + AuqMI . (8) 



In the limit d ^ 3 we then obtain the singular part of the bare Helmholtz free 
energy as a finite function r(rQ,Mo,Mo) whose two-loop expression has been 
derived recently (see Eq. (31) of Ref. [6]). 

We have been able to perform an analytic calculation of the three-loop contri- 
butions to r(rQ,-uo, Mq) near the coexistence curve below and for Mg = 
above T^. The details of this calculation will be presented elsewhere. Some 
comments are given in Appendix A. The resulting bare perturbation expres- 
sion for the Helmholtz free energy at Tq 7^ reads up to three-loop order at 
d = 3 



f(r^„«o,Mo) 



1 / 



3 b-l 1 



T,Ml + uoM', + E E E (-1)' 2-'-'=F,^fc(«;, n) 



6=1 1=0 k=0 

4- 



-21 



(24mo)^ 



In 



(24mo)^ 



(9) 
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where 



w{ro,Uo,Mo) = ror/roL 



(10) 



is a non-perturbative parameter. Above Tc its largest possible value is w{rQ, uq, 0) 
1 (when Tq > and Mq = corresponding to T > Tc at Hq = 0). The small- 
est possible value of w is attained on the coexistence curve well below T^. 
where, at given Tq < and to lowest order in uq, w{rQ,uo, Mo{rQ,uo)) — 
37r~^-uo(— 2rQ)~-'^/^ + O('Uq) as follows from Eq. (22) in Sect. 3. The coefficients 
Fbikiw, n) depend on w and n and determine the contributions in 6-loop order. 
The one- and two-loop coefficients are [6] 



Fioo(ty,n) = 



1 

3847r2 
1 



1 + (n - l)w^/^' 

3 + 2(n - l)w^/2 + (n^ - l)w 

, l + 2^x;V2 
n — 1) in , 



2887r2 

The new three-loop coefficients F300 and F301 read in analytic form 



(11) 
(12) 

(13) 
(14) 



F3oo(w,n) 



1 f 3 

ns + 24 In - - (n - 1) 

184 327r3 1 4 ^ ^ 



-1/2 + 2n - 6 



+ 8 In h w^'^ - 6n - 9 + 4(n + 1) In -— 



+ 81n ^ J +it;(n-l) 

^ 23^ {3 + (n - 1) [1 + (n + 2)^^/^] } 



(15) 
(16) 



In calculating the coefficients F310 and F320 below Tc we have confined ourselves 
to the vicinity of the coexistence curve where an expansion with respect to w 
is justified. For the purpose of calculating the order parameter this expansion 
must include the terms of linear order in w whereas for the calculation of the 
specific heat all terms can be neglected that vanish in the limit w 0. The 
latter statement applies to T300 and T301 as well. The result for the coefficients 
F310 and F320 is 



310 



276i8^{^"'-^^ + ^°^^^^(-0-^"-^) 



4iD-^/2 
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+ 4n + 2 - (n + 2)7r^ - 12 Lia (|) - 32 In 2 - 6 (In 3)^ 
+ ty^/^ (lOn + 32 - 16 (2n + 3) In 2 + 48 In 3 - 8(n + 1) In 

+ -^^(8471 - 100 - 1281n2) + 0(w^/^ -ly^/Mn w) 



w 



(17) 



F32o{w,n) 



|4321n 1 - 324Li2(-i) - 432ci - 27n' 



1658887r3 
-in-l 



16^-1/2 ^ 3ri + 14^2 ^ i8(in3)2 + 36 Lis (|) 



+ 16(c2 + 4Li2(-2) -2Li2(-i) + 



61n3-ln2 - — 
3 



In 2 



^ + 16iD^/^(7 - n + (n + 1) In (16w) + 2 In 2 - 6 In 3) 



224 



+ 4w |^4c2 - 12n ^ 



+ 7r^ + 6 



61n3-ln2- 



16 
15 



In 2 



+ 12 



2Li2(-2) - Li2(-i)] ]+0{w^/\w^/Hiiw) 



(18) 



Here Li2(a;) = J^t ^ In (1 — t) dt is the dilogarithmic function. The constants 
Cl and C2 are given by 



± 

I 





dx 



VQ - 2a;2 



,3 ,3 + a; x / 3 + a; x ^ 2 + x 
m - + in 1 m 1 in 



4 2 + a;2 + xV 3 2-x 4 



0.0217376333 



(19) 



and 



C2 



dx 



AV2 ' '"^ y/TTx- 



, X ln(l + a;) 

In \ ^ - 

1 + X X 



= 0.973771427. (20) 



Above Tc, at Mq = 0, the coefficients in Eqs. (11)-(16) are taken at w = 1 
and the coefficients F310 and F320 do not contribute. All other coefficients 
F'iik vanish. Eqs. (9)-(20) are the main result of our paper. They provide the 
basis for deriving the analytic form of the amplitude functions of the order 
parameter and the specific heat for general n. For the case n = 1, Eqs. (9)- 
(20) reduce to the three-loop part of Eq. (3.18) of Bagnuls et al. [7] and of 
Eq. (3.3) of Halfkann and Dohm [8] and agree with the numerical values of 
their coefficients Ffji^ for 6 = 1, 6 = 2 and 6 = 3. The latter agree also with 
those in Eq. (Al.l) of Guida and Zinn- Justin [9] for n = 1. 

Eq. (9) contains logarithms of the coupling uq as expected because of the non- 
analytic uq dependence of (5ro, Eq. (6). Furthermore, the terms proportional 
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to n — 1 depend non-analytically on Tot through w~^l'^ ^ w^l'^ and w^/^lnw 
(and higher orders in the neglected terms). These non-analyticities will lead 
to perturbative terms in the derivatives of the free energy (with respect to 
Mo) that diverge when the coexistence curve is approached (T < Tc, /to — >^ 0). 
This is again the effect of the Goldstone modes that was found previously in 
two- loop order [6]. For 0(n) symmetric quantities, however, these divergences 
should cancel among themselves [34]. This is indeed the case, at least up to 
three-loop order, for the square of the order parameter, the Gibbs free energy 
and for the specific heat as we shall see below. 



3 Bcire order pcirameter 

The order parameter Mo(rQ, UQ,ho) is determined by inverting the equation of 
state, 

ho{r'^, uo, Mo) = (r(„ Uo, Mo) . (21) 



This inversion should be performed itcratively at finite ho ^ 0. Several per- 
turbative terms of dt / OMq exhibit Goldstone singularities for /iq — > below 
Tc which arise from the three-loop terms ~ w"^/'^, w^/'^ and w^l'^\n.w (after 
an expansion of w dependent terms with respect to Uq at /lo 7^ 0), in addition 
to the known [6] two-loop singularities, as well as from the expansion of one- 
and two-loop terms (with respect to Uq contained in M^{r'f^,U(), ho) at finite 
ho). We have verified that all Goldstone divergences cancel among themselves 
which leads to a finite result for Mq as a function of Tq and uq in the limit 
ho ^ at d — 3. This function reads for Tq < 



+ 



+ 



10-n + 4(n- l)ln3-2(n + 2)ln■ 



_Uo_ 
87r2 



-2r' 



(24iio)2 

2736n - 5904 - 6480ci + 240(n - l)c2 
(n-l)Li2(|)+9Li2(-|); 



19207r3 [ 
- (7571^ - 5n + 875)7r2 - 1260 
+ 960(n- 1) [2Li2(-2) - Li2(-|)] - 630(n - l)(ln 3) 



10(n- l)ln2-60(n- l)ln3 + llln-561 

2r' 1 



-481n2 

+ 240(12n - 57) In 3 - 1440(n + 2) In 



(24xio)2 



9 



+ O{ul,ullauo) 



(22) 



It contains the expected [6] logarithmic Uq dependence ~ ullnuo at three- 
loop order. The logarithmic terms in uq can be absorbed by employing the 
pseudo-correlation length [5,6] ^_ as a temperature variable below instead 
of Tq < 0. The relation between Tq and ^_ is determined by the two-point 
vertex function and is given up to three-loop order by (see Appendix B) 



- 2r' = C 



-2j 1 , "- + 2 ^ n + 2 2 



+ 



TT 

n + 2 
IOSttS 



1385 
~W8 



+ 41n(24Moe-) 



1314n + 13047 + 576(n + 8) Lig (-| 



+ 48(n + 8)7r^ + 8 (43n + 182) In - 



+ OiutC) . (23) 



This leads to the bare perturbative expression for the square of the order pa- 
rameter Mq (^_, Uq, 3) that is finite in three dimensions and is free of logarithms 

of Uq, 



M2(e_, «o, 3) = CM E a^m{M-r-' + o{uie) 

lm=0 



(24) 



with the coefficients up to two- loop order [6], 



a 



_ 1 
<.o-g, 

1 



^ n-l)ln3--l-,(1169n + 1042), 



27^2^'" 8647r2 
and the new three-loop coefficient 



(25) 
(26) 
(27) 



= -.^oL J 36 (565n2 + 5056n + 7744) - 58320ci + 2160(n - l)c2 
1728071'^ [ 

+ 157r^(19n^ + 643n + 499) + 180 (64n^ + 640n + 457) Li2(-|) 

- 11340(n- l)Li2(|) +8640(n- 1) [2Li2(-2) - Li2(-i)" 

+ 80 (86n^ + 860n - 811) In 3 - 16 (860^^ + 8357n - 7867) In 2 

- 5670(n - 1) (In 3)^ + 4320(n - 1) (6 In 3 - In 2) In 2 i . (28) 
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For n = 1, Eqs. (24)-(28) agree with Eq. (3.13) and the numerical values in 
Table 2 of Halfkann and Dohm [8]. 



4 Bare Gibbs free energy 



The bare Gibbs free energy f is determined by the bare Helmholtz free energy 
and the order parameter as 

.F(ro, Wo, /io) = r(ro, uq, Mo(ro, wq, h^)) - hoMo{rQ, uq, ho) . (29) 



Above Tc we obtain the Gibbs free energy ^+(ro, Uq) at ho — by substituting 
Mo = and = 1 into f . Prom Eqs. (9)-(16) we find at d = 3 



:F+(ro,uo)^r(ro,uo,0) 



n /3/2 n{n + 2) , 2n{n + 2) ^ ,1/2 
^0 rA-\-) ""0^0 ~ 



127r 



2 '"o'^o 



81n^ + 41n-J^- 



(47r)2 ^"'^ (47r) 

+ 0{ul,ullnuo) 



n — 6 



(30) 



for Tq > 0. Below we make a perturbative expansion of J^, Eq. (29), in uq 
at /iq 7^ and take the limit 



(ro, -Uq) = lim :r(ro, Mq, /lo) , < . 



(31) 



The result is at d = 3 



{-2r',f {-2r'o)'/' 



6Auo 127r 

Mo(-2rd 



+ 



167r2 

«g(-2rd)V2 



6 + 2(n- l)ln3- (n + 2) In ■ 



-2r' 



3847r3 



(24ko)2 

16 (lln + 7) - 1296ci - 48(n - l)c2 

+ (971^ + n + 17)tt^ + 36(n - 1) Li2(|) + 324 Li2(-|) 
+ 96(n- 1) [Li2(-|) -2Li2(-2)] +18(n- l)(ln3)' 
+ 48(n- l)(ln2-61n3)ln2-48 (4n + 17) ln3 

2r' 



+ 16 (31n + 95) In 2 + 96(n + 2) In 
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+ O{ul,ullnuo) . 



(32) 



As expected, all Goldstone divergences cancel among themselves in this limit. 
The terms of 0{w^/^), ©(w^/^lnw) and 0{w) in Eqs. (15)-(18) do not con- 
tribute to the free energy ^-(rg, Uq) and to the specific heat (9^JF_ (rg, uq) / {dr^Y 
on the coexistence curve because these terms vanish in the limit /iq — ^ 
(iD — > 0). Therefore, for the purpose of an application to the specific heat in 
Ref. [10], it was sufficient to evaluate most of the diagrams in Fig. 1 directly 
at Tor = (iD = 0), see also App. A. This simphfication, however, is not 
applicable to the calculation of the order parameter performed in the present 
paper. 

We note that the function P-{rQ, Uq) can also be considered as the Helmholtz 
free energy r(rQ, Uq, Mq) below T^. in the limit where Mq approaches the spon- 
taneous value Mo(rQ,Mo), Eq. (22), of the order parameter on the coexistence 
curve. Thus the Helmholtz free energy remains finite in this limit, as expected 
[34] . In order to express J^± in terms of the correlation lengths ^± we use (see 
Appendix B) 



, ^-2 f -, n + 2 . n + 2 . ^ 



+ 



TT 

n + 2 



27 



+ 21n(24iioC+) 



3 (3n + 22) - 144(n + 8) Li2(-|) 
3" 



547r3 

12(n + 8)71^ - 2 (43n + 182) In 



(33) 



for Tq > and ro(C_,Mo) as given in Eq. (23) for Tq < 0. (We recaU that, for 
d — 3, the right-hand side of Eq. (33) contains a logarithmic Uq dependence 
only at 0{ul), as can be seen in Eq. (4.5) of Ref. [20] for e = 1. The same 
comment applies to Eq. (23).) Substituting Eqs. (33) and (22) into Eqs. (30) 
and (32) yields the Gibbs free energy up to three-loop order 



-3 . 



n 

~ 12^ ' 

+ \n — 

9 . 



3n{n + 2) n{n + 2) 

(«oe+) - 

{uo^+f + O{ulil,ul^l\n{uoi 




(34) 



and 



(35) 



with the coefficients up to two-loop order 



12 



a_Q — 



a_2 



1 

64' 
1 

"96n 
1 



34 567r2 

and the new three-loop coefficient 



(3n + 14) , 

54n2 - 737n - 394 + 432(n - 1) In 3 



(36) 
(37) 
(38) 



1 



17971^ - 3875n - 10626 + (33^^ - 471n - 615)7r^ 



34 567r3 

- 11664ci - 432(n - l)c2 - 12 (48n^ + 480n + 525) Li2(-|) 
+ 540(n - 1) Li2(-i) - 1728(n - 1) Li2(-2) 



+ 54(n - 1) 



51n2-42 ln3 



In 2 + 16 (43n2 + 547n + 1219) In 2 



-8(97n^ + 538 + 1174)ln3 



(39) 



For n = 1 this Gibbs free energy is given for T < Tc in Eqs. (3.14), (3.15) 
and Table 2 of Ref. [8] where it is denoted by r^o{^^,Uo). The logarithmic 
contributions 0{ullii{uo^±)) in Eqs. (34) and (35) correspond to Eq. (3.15) of 
Ref. [8] and are due to a specific d — 3 divergent diagram arising only at four- 
loop order [7] . Its pole term is to be subtracted within the overall subtractions 
mentioned in Section 2. 



5 Bare specific heat 

The specific heat abve (-I-) and below (— ) Tc at ho — is determined by [5,6] 

C^^CB-alt^l''\r'„Uo) (40) 



where 



± 



(41) 



Cb is a. noncritical background value and Gq is a constant defined by ro — Tqc 
aot where t ^ {T - T^)/T^. The result is 



n /-1/2 , n(n + 2) ^ ,.3/2 



167r 



2(47r)3 ""^'^ 



n — 6 — 8 In 
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+ 41n 



+ 0(M^,'u^lnMo) 



(42) 



and 



4(-2r(,)-3/2 



3847r3 



16 (lln + 7) - 1296ci - 48(n - l)c2 

+ (9n' + n + 17)7r' + 36(n - 1) Lis (|) + 324 Lis (-|) 
+ 96(n-l) [Li2(-|) -2Li2(-2)] + 18(n - l)(ln3)' 
+ 48(n- l)(ln2 -61n3)ln2 -48 (4n+ 17)ln3 

— 2r' 

+ 16 (31n + 95) In 2 + 96(n + 2) In ■ 



(24«o)^ 



+ 0{ul,ulhiuo) . 



(43) 



Unlike the corresponding two- loop result [6] for f±'°\ the three-loop terms 
of f 1^'°'' contain a logarithmic -UQ-dcpcndcncc (as long as they are considered 
as a function of Tq). This can be absorbed by turning to the formulation in 
terms of the correlation lengths ^± [5,20]. The relation between Tq and ^± is 
given in Eqs. (33) and (23). This leads to the bare perturbative expression 
for the specific heat at /lo = that is finite in three dimensions and is free of 
logarithms in uq, 



n n{n + 2) 



167r 327r2 



27n + 160 + 1081n- 
4. 



+ 0{ulil) 



(44) 



and 



fL^'°)(e-,x.o,3) \ E a'Z\uo£,-r-' + 0{uli 

lm=0 



3f3 ' 



(45) 



with the coefficients up to two- loop order [6] , 



,(2,0) 



„(2,0) 



1 

'4^' 



(46) 
(47) 



14 



(48) 



and the new three-loop coefficient 



.(2,0) 



34 567r3 



1188n^ + 11876n + 16840 - 11664ci - 432(n - l)c2 



+ 97r^(9n^ + n + 17) + 324(n - 1) Li2(|) + 2916Li2(-|) 
+ 864(n- 1) [Li2(-i) - 2Li2(-2)] + 162(n - l)(ln3)2 
+ 432(n- l)(ln2 - 61n3) ln2 -432 (4n+ 17)ln3 

+ 144 (31n + 95)ln2 



(49) 



For the case n = 1, Eqs. (45)-(49) agree with Eq. (3.17) and the numerical 
values in Table 2 of Halfkann and Dohm [8] . 



6 Renormalization and amplitude functions in three dimensions 



In the preceding Sections no renormalizations have been used, except for the 
shift of the temperature variable Tq by Stq, Eq. (6), and for the additive sub- 
traction of the Helmholtz free energy. All other quantities are unrenormalized. 
Because of the super- renormalizabihty of the 0^ theory for d < 4, the shift of 
Stq and the additive subtraction are sufficient to make the order parameter 
Mq{C,-,Uq,3) and the vertex function r±.''^\$^±,UQ,3) finite at infinite cutoff 
in three dimensions as long as ^± is finite. The resulting bare perturbation 
series in Eqs. (24), (44) and (45), however, are obviously not applicable near 
Tc where uq^± diverges, thus the series have to be mapped from the critical 
to the non-critical region. 

This mapping is achieved by turning to the renormalized theory as defined 
below and by introducing the renormalization scale /i which can be varied 
via the renormalization-group equation (RGE) [35]. The solution of the RGE 
implies a decomposition of thermodynamic quantities into a product of am- 
plitude functions and exponential parts. This decomposition is most natural 
and simple within the minimal subtraction scheme [36] where the exponen- 
tial parts are determined entirely from the pure d — A pole terms ~ e~" of 
the renormalization constants Zr{u, e), Zu{u, e), Z^{u, e) and A{u, e) [5,20]. As 
shown previously [5,17,20], the use of these pure pole terms does not imply the 
necessity of working near four dimensions and of using the e-expansion. This 
is a nontrivial aspect of the minimal subtraction scheme that is not widely 
appreciated in the field-theoretic literature. Here we make use of the fact that 
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the minimal subtraction scheme can well be combined with the concept of a 
(f)^ theory at fixed dimension d < 4. 

Our d — 3 approach differs from the d — 3 approach of previous work [7,9,21- 
24] where renormalization conditions were used. These renormalization con- 
ditions are different above and below Tc whereas the minimal renormalization 
applies to both T > Tc and T < Tc in the same form. Within our scheme the 
minimally renormalized quantities are introduced in d dimensions as 



r = ^(ro - roc) , u = fi ^AaZ^ ^Z^uo , $= ^^Vo , (50) 
M\i_,u,ii,d) = Z^'Ml{i_,ii^A-/Z^Z^^u,d) , (51) 

Tf'\i±,u,ii,d) ^ Zltt^'\i±,lx^Z^Z^^Atu,d) - -^ii-^AaA{u,e), (52) 
where ^± are the correlation lengths defined in Refs. [5,20] and 

r(3-d/2) 



is a convenient geometric factor [17]. We recall that the amplitude functions 
to be defined in Eqs. (54) and (55) depend on the choice of this geometric 
factor. In our calculations, Eqs. (50)-(53) are used directly at d = 3. The 
analytic form of the renormalization constants Zr{u, e), Zu{u, e), Z<^(m, e) and 
A{u,e) is given in Eqs. (2.13), (2.16)-(2.19) and (B1)-(B18) of Ref. [10] for 
general n up to five-loop order. As is well known, the minimal renormaliza- 
tion constants in Eqs. (50)-(52) have the property of absorbing the remaining 
ultraviolet divergences for c? — >^ 4 that were not absorbed by the mass shift 
and subtractions mentioned in Sect. 2. In the context of our d — 3 theory this 
property is, of course, irrelevant. The crucial aspect of the renormalization of 
Eqs. (50)-(52) at e = 1 is that they provide the mapping and decomposition 
mentioned above via the integration of the RGE for the renormalized quanti- 
ties and T^^'^\ The critical behavior of these quantities then evolves from 
the infrared divergences of the Z-factors Zi{u, 1) as u ^ u* [20,23] and from 
the singular behavior of A{u, 1) as it — > u* [26]. 

Dimensionless amplitude functions and F_ of the renormalized order 

parameter and specific heat can be defined at fixed dimension 2 < d < 4 
according to 

U{^^i-,u, d) - e-^M\i_, u, /X, d) (54) 

and 

F±{f^^±, u, d) = -4/xM^iri''°)(e±, u, /X, d) . (55) 
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These functions remain finite also in the limit 4 (at finite ^-t) [5,20] 

provided that the appropriate subtractions (including all additive pole terms 
for (i < 4) have been performed in Sect. 2. In the application of the solution of 
the RGE the amphtude functions appear in the form [5] with the non-critical 
arguments ii^± — 1 in three dimensions 

U{u) = U{l,u,3), F±(«) = F±(1,«,3). (56) 



These functions should be smooth and well behaved in the entire region < 
u < u* [20] . They have the power series 



Mu) = i £ c-^«™ , (57) 

" m=0 

oo 

m=0 
-1 oo 

" m=0 

with n-dependent coefficients. These series have a zero radius of convergence 
but are (presumably) Borel resummable [20,28] (see Sect. 7). Note that the co- 
efficients of these series depend on the choice of the geometric factor A^^. Our 
choice, Eq. (53), minimizes the explicit dimensional dependence of the lowest 
order coefficients in Eqs. (57)-(59) [5,17,20,37] and is expected to improve the 
convergence properties of the series. This is of relevance in the context of Borel 
resummations based on low-order information. 

From Eqs. (24), (44) and (45) we obtain the analytic expressions of the coef- 
ficients for the order parameter 



1 

c^i = 0, (61) 



c.-o-^, (60) 



= 7^(160 - 82n) + -(n - 1) In 3 , (62) 



^(160- 82n) + -i 

:-o = 2500n2 + 65104n + 29056 + 8640 (5n + 22)C(3) + 58320ci 

10807r[ ^ ^ 

+ 2160(n - 1) [4Li2(-|) - ca] - 157r^(19n^ + 643n + 499) 

- 180 (64n^ + 640n + 457) Li2(-|) + 11340(n - 1) Li2(|) 

- 17280(n - 1) Li2(-2) + 5670(n - l)(ln 3)^ + 4320(n - l){\n2f 

- 25920(n - 1) (In 2) (In 3) - 80 (194n^ + 1616n - 1675) In 3 
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+ 16 (860/1^ + 8357n - 7867) In 2 | , (63) 
for the specific heat above Tc, 

4o = > (64) 

= -2n(n + 2) , (65) 



7 .,4 

n h 4 m - 

27 3 



= -4n(n + 2) 
and for the specific heat below Tc, 



(66) 



Cfo = ^ > (67) 

c^i = -4, (68) 
c-2 = 8(10-n) , (69) 

c-3 = --I(i080n2 + 3464n + 31120) - 128 (5n + 22)C(3) - 864ci 

2 

- 32(n - l)c2 + y (1871^ + 2n + 34) + 24(n - 1) Li2(|) 

+ 216Li2(-i) + 64(n - 1) Li2(-|) - 128(n - 1) Lis (-2) 

+ 12(n - 1) (In 3)2 + 32(n - l)(ln2)2 - 192(n - 1) (In 2) (In 3) 

32 

-32(4n + 17)ln3 + — (31n + 95)ln2. (70) 

3 

The terms up to 0{u) are the previous two-loop results [5,6]. In Table 1 we 
give the numerical values of the coefficients c~^, c^^ and Cp^ up to three- loop 
order for n = 1, 2, 3 as obtained from Eqs. (60)-(70). 

For n — 1 the numerical values of the coefficients given previously [8] 
agree with ours up to eight digits. The difference in the ninth digit is due to 
a numerical inaccuracy of the Z-factors in three dimensions that were avail- 
able previously only in numerical form [38]. (A similar comments applies to 
the numerical values of the series for n = 1 [8] mentioned in the preceding 
Sections.) In the present work we have used the Z-factors in the analytic form 
for general n as given in Ref. [10]. 

Numerical results of a six-loop calculation (above for n = 1, 2, 3) and of 
a five-loop calculation (below Tc for n = 1) of the diagrams contributing to 
F± in three dimensions have been presented by Baker et al. [22], by Bagnuls 
and Bervillier [23] and by Bagnuls et al. [7], respectively. On this basis the 
coefficients Cp^ for n = 1, 2, 3 and Cp„^ for n = 1 were calculated in numerical 
form by Krause et al. [38] and by Halfkann and Dohm [8] who used a two- 
loop approximation for the additive renormalization constant A{u,e). The 
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Table 1 

Coefficients of f^{u) according to Eqs. (57), (60)-(63), and Cp^ of and 
F-[u) according to Eqs. (58), (64)-(66) and (59), (67)-(70), respectively, for n = 1, 
2, 3 up to three-loop order. For and c^^, m refers to u'^~^ corresponding to 
m-loop order whereas for c^^, m refers to corresponding to {m+ l)-loop order. 







m 








n = 


1 





(327r)-i 


-1 


0.5 






1 





-6 


-4 






2 


0.919561893 


-22.6976284 


72 






3 


-80.8015258 




-5189.75474 


n = 


2 





(327r)-i 


-2 


0.5 






1 





-16 


-4 






2 


0.652241285 


-92.5270090 


64 






3 


-83.0064428 




-5918.07320 


n = 


3 





(327r)-i 


-3 


0.5 






1 





-30 


-4 






2 


0.384920676 


-233.488142 


56 






3 


-82.6969869 




-6607.95641 



corrected coefficients for Cp^ ior n — 1, 2, 3 and Cp^ for n = 1 using A{u,e) 
in five-loop order have been presented recently in numerical form [10]. Our 
present results, Eqs. (64)-(70), provide the analytic form of these coefficients 
up to three-loop order for general n. The three-loop coefficient Cp2 has been 
calculated independently by Burnett [39]. 



7 Results and discussion 

Within the minimally renormalized 0^ field theory in three dimensions we 
have derived the three-loop contributions to the amplitude functions /</,(«) 
of the square of the order parameter, Eqs. (57), (63), and to the amplitude 
functions F+{u) and F_{u) of the specific heat above and below T^., Eqs. (58), 
(66) and (59), (70) for general n. Similar to the previous results below 
for n = 1 [7,8,10], the coefficients c^^ and Cp^ for n = 2 and n = 3 given in 
Table 1 have alternating signs and increase considerably in magnitude. Clearly 
a resummation of these series is necessary in order to obtain quantitatively 
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reliable results. For a description of the method of Borel resummation we refer 
to Refs. [20,28] and, in the present context, to Ref. [10]. 



While the previous two-loop results [6] did not yet provide sufficient infor- 
mation for a controlled resummation procedure and thus did not yet lead to 
an error estimate, the new information on c~3 and c^^ presented here makes 
possible to perform the Borel resummation for n > 1 below Tc with reasonably 
small error bars. The reliability of this Borel resummation based on three-loop 
results (with four low-order coefficients, see Table 1) has been demonstrated 
in detail in a related paper [10]. We refer to this paper for a description of the 
method of determining the error bars. For the case n = 2 (corresponding to 
the superfluid transition of "^He) the resummation result is 



at the fixed point m* = 0.0362 [10]. In addition, using the new coefficient 
iov n — 2 and n = 3, a Borel resummation of the fixed point value of the 
amphtude function f^{u*) has been performed by Monnigmann [40]. (For this 
case the Borel resummation parameters h and a [20,28] are h = 3.5+n/2 = 4.5 
and -0.43 < a < 0.47 for n = 2, and 6 = 5, -0.62 < a < 0.71 for n = 3.) 
The result is 



at the fixed point [40]. Thus the deviation from the zero- loop term c^q — 
(327r)~^ 0.009947 is very small, confirming the previous expectation [5]. 

The amplitude functions are of physical relevance not only at the fixed point 
but also in a finite interval < u < u* corresponding to a finite distance 
from criticality where non-asymptotic corrections become non- negligible [18]. 
Since the Borel resummed amplitude functions F-{u) and /</,(«) are smooth 
and slowly varying functions for n = 1 (see Figs. 1 and 4 of Ref. [8]) the same 
property should hold also forn = 2 and 3. Empirical evidence in support of this 
fact comes from analyses of experimental data of ^He in the non-asymptotic 
region [16,18,41]. Therefore, within the accuracy of the present results, these 
functions can be approximated by the simple representations [8,38] 



u*F^{u*) = 0.384 ± 0.025 



(71) 




n = 2, 
n = 3 



(72) 
(73) 



F_{u) = {2u) 



-1 



4:{l+dFU) 



(74) 



and 



(75) 
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with A3 = (47r)~^. Here the lowest-order terms are chosen to coincide with 
the one-loop results. The coefficients dp and ci^ are determined by requiring 
that Eqs. (74) and (75) agree with the Borel-resummed fixed point values of 
Eqs. (71) and (72). This implies 

dp = -5.49 , d^ = 0.422 (76) 



for n = 2. In view of the present error bars, a more refined representation 
of F^{u) and /</,(«) does not seem to be warranted until four-loop (Borel- 
resummed) results with smaller error bars become available. 



0.008 



1 ' 1 ' 1 ' 

2 loop 

1 






& 1 loop 




n=2 






Borel 






1 


3 loop 

1 



0.00 0.02 0.04 

u 

Fig. 2. Partial sums fjp^\u) = \^Y^m^{) (^m'^^™' of the amplitude function 
f(f,{u) = /0(l,tt, 3), Eq. (57), for the square of the order parameter in three di- 
mensions for n = 2 multiplied by it, as a function of the renormalized coupling u 
from M = (zero- loop order) to M = 3 (three-loop order) (dashed lines). The 
solid curve is Eq. (75) with = 0.422 representing the Borel summation result. 
The error bar indicates the error of the Borel resummed fixed point value /^(tt*) at 
u* = 0.0362, Eq. (72). 

In Figs. 2 and 3 we have plotted /^(m) and F_{u) for n = 2 according to the 
representations in Eqs. (74) and (75) (solid lines), together with the zero-, one- 
, two- and three- loop approximations in the range < it < it* (dashed lines). 
The amplitude functions f^{u) and F±{u) enable us to calculate and predict 
(rather than fit) the various ratios of amplitudes of leading and subleading 
terms of the critical temperature dependence of the order parameter and the 
specific heat [5]. These ratios are universal quantities [13]. 

As a first application we consider the specific heat which plays an important 
role in testing the renormalization-group prediction of critical-point universal- 
ity along the A-line of ^He [11,12]. Combining the three-loop result for F_{u*), 
Eq. (71), with the recent [10] five- loop result for the additive renormalization 
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0.6 



uFJu) 



0.4 



0.2 



0.0 




n=2 



Borel 



0.02 



1 loop 



3 loop 



0.04 



Fig. 3. Partial sums F^\u) = ^ Z^m=o <^Fm^"^ amplitude function 

F^{u) = 3), Eq. (59), for the specific heat below in three dimensions 

for n = 2 multiplied by function of the renormalized coupling u from M = 

(zero-loop order) to M = 3 (three-loop order) (dashed lines). The solid curve is 
Eq. (74) with dp = —5.49 representing the Borel summation result. The error bar 
indicates the error of the Borel resummed fixed point value F_(u*) at u* = 0.0362, 
Eq. (71). 

constant A{u, e) leads to the universal amplitude ratio [10] 



A+ 

— = 1.056 ±0.004 

A- 



(77) 



for the asymptotic specific heat in three dimensions. This is in excellent agree- 
ment with the high-precision experimental result A'^/A~ — 1.054 ± 0.001 for 
^He near the superfiuid transition obtained from a recent experiment in space 
[15]. Our result, Eq. (77), is a significant improvement over the previous pre- 
diction [4] A+/A- = 1.0294 ± 0.0134 obtained from the e = 4 - d expansion 
up to 0{e^) [42]. 

As a second application we employ our Borel-rcsummed three-loop result for 
to calculate the universal combination of amplitudes [13] 

^~ Al, (3/2 - 2i/P^)2/3 ■ IGtt y^'t'^ J ^ ' 



where and Am are the leading amplitudes of the susceptibility above and 
of the order parameter below Tc, respectively [13]. The structure of the right- 
handed side of Eq. (78) follows from previous work [5,38]. Here = /</,(«*) is 
determined by Eqs. (72) and (73). The quantities P* = P(m*), F* = F+(m*), 
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B* = B{u*) and = f(^\u*) have been calculated in Refs. [10] and [38], 
respectively. The values for the critical exponents a and u are given in Ref. [10], 
the value of /3 is taken from Ref. [43]. The resulting values for Rc in three 
dimensions are 



These results are significantly more accurate than the e-expansion values [29] 
up to O(e^) whose extrapolation to e = 1 suffers from considerable ambiguities 
[29]. The numerical estimate Rc = 0.165 of high-temperature series expansions 
for n = 3 [44] is not far from our value, Eq. (80). Additional numerical and 
experimental studies would be desirable. 

Further applications of our results will be given elsewhere. 
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A Three-loop vacuum diagrams 

In this Appendix we present some details of the calculation (at infinite cutoff) 
of the three-loop vacuum diagrams shown in Fig. 1 in three dimensions. There 
are six different topologies of vacuum diagrams (denoted as (A), (B), . . . (F) 
in Ref. [27]) whose integral expressions are given by 



i?c = 0.123 ± 0.003, n = 2, 
i?c = 0.189 ± 0.009, n = 3. 



(79) 
(80) 




P1P2P3 



(m^+p^)(m^ + lPi-P2|2)(m^+p^)(m^+|p2-p3|2)(m^+p^)(m^ + lP3-pJ2) , 

(A.l) 



1 



1 2 3 4 = 




P1P2P3 



/// 




/// 



(m^+p^)(m^+|P3-PiP)(m^+|P3-P2P)(m^+p^)(m^+p^) 



1 



(A.3) 



P1P2P3 
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2 



(m2+pf)(m|+pi)(m| + |p,+P2|2)(m|+|pi+P2|2)(m|+p§) 



|2Vr„.2_i_„2-| , 



P1P2P3 




(mJ+p^)(m^+p^)(m^+p^)(m^+|Pi+P2+P3P) 



21 ' 



1 



ayV_y 7 7 7 ("I? +Pl ) ("*2 +P2 ) ("*3 +P2 ) ("*4 +^3 ) 

P1P2P3 



(A.4) 
(A.5) 
(A.6) 



In our case, the masses mf are given by tql or tqt of Eq. (8). Three of the dia- 
grams in Fig. 1 arc infrared divergent in the hmit Tot ^ in three dimensions: 



type (B) with mf = mi = Tql and m 



^otj type (D) with 



rrin 



Tql and m 



Tqt and type (F) with ml 



roL 

and m| = 777.3 ~ '"ot- The leading tot dependence of these diagrams appears 
in the terms of 0{w~^^'^) in Eqs. (15), (17) and (18). All other diagrams in 
Fig. 1 are free of infrared divergences in the limit tot ^ at d = 3. 

Most difficult is the evaluation of the "Mercedes" diagram (A), Eq. (A.l). It 
has been calculated by Rajantie [27] in three dimensions only for two special 
cases (i) where one of the masses mf vanishes, or (ii) where all the masses are 
equal. The case (ii) corresponds to our first Mercedes diagram in Fig. 1 with 
= f^oL and is given by Eq. (55) of Ref. [27]. The last two Mercedes diagrams 
in Fig. 1, however, contain the two different (longitudinal and transverse) 
masses Tql and Tot- This case is considerably more difficult than the cases (i) 
and (ii) mentioned above. We have been able to calculate these diagrams for 
small w — tqt/tql > in three dimensions as 




-3/2 
(47r)3 



^Li,(-2) + ^(ln3)(ln2) + ^ + 3^^/2ln^ 
+ |3Li2(-2) + 3(ln3)(ln2) +41n2 + ^ 



+ C>(#/2,w3/2lni(;) 



-3/2 



OL 



C2 + 2w^/'^{\nw + 61n2 - 2) + w{c2 - 6 - 41n2) 



(A.7) 



(A.8) 



where C2 given by Eq. (20). 

For finite mf the diagrams (B)-(D) and (F) are finite in three dimensions at 
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infinite cutoff whereas (E), Eq. (A. 5), has a single pole ~ (rf — 3) ^. Analytic 
results for the diagrams (B)-(D) and (F), Eqs. (A.2)-(A.4) and (A. 6), are 
given in Eqs. (A.24), (A.25) and (A.27) of Ref. [27] and in Eq. (5) of the 
Erratum [27]. 

Regarding the d = 3 divergent diagram (E), the MS scheme was used in 
Ref. [27] which differs from our approach. Within our approach the result for 
diagram (E) near d = 3 is represented as 




u 



Uq i( 2{ 3\ 4) 



-1+d/e 4 

° E 

i=l 



rrii 



— '- 2 + -7 - - n 

e- 1 4 ^ 4 



+ iln^ + ln^ + 0(e-l) 

2 Uq Uq 



(A.9) 



with 7 = CeuIci- being Euler's constant and e = 4 — (i. The pole term ~ (e — !)"''■ 
is cancelled by a term of order 0{u^) that arises from performing the mass 
shift ro = Tg + 5rQ in the one-loop contribution of the Helmholtz free energy. 
Eq. (A.9) corresponds to Eqs. (21) and (A.26) of Ref. [27] and Eq. (1) of the 
Erratum [27]. Note that in our Eq. (A.9) there is no additional parameter 
corresponding to the MS scale parameter // of Ref. [27] . 



B Correlation lengths in three-loop order 



In this Appendix we derive Eqs. (33) and (23). Above Tc, the correlation length 
is defined via 



q=0 



(B.l) 



where x+{q) ^ = T'Q''^\q, Tq, Uq) is the inverse susceptibility at finite wavenum- 
ber q. The two-point vertex function ro°'^''(g, Tq, Mq) is given by 



ro°'^Hg,ro,'Uo) = ro + - Eo(g,ro,Mo) 
where the self-energy 

oo 

Eo(g,ro,'Uo) = J2i-'^oT^o^\q,ro) 



(B.2) 



(B.3) 



771=1 



is the sum of all one-particle irreducible m-loop diagrams with two (ampu- 
tated) external legs. In Eqs. (A.1)-(A.5) of Ref. [6] the functions r'Q''^\q, tq, uq), 
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Eq. (B.2), and ^+(ro, Uq), Eq. (B.l), arc given in their two- loop approximation. 
The diagrams of the three-loop contribution to Eq. (B.3) are given by 



+ 384(n + 2) 




+ 64(n + 2)^ 



+ 128(n + 2)^ 




+ 128(n + 2){n + 8) .(B.4) 



For g = the lines denote the standard propagator (ro + ^ above Tc. The 
derivative of Eq^'' (q, ro) with respect to yields 



q=0 



384(n + 2f 



e d 



d dml 
128(n + 2)(n + 8) 




2ro 92 

+ 



e d rK\ 2ro 

U 2/ 3\ 41 -|- 
_5, 



m|=ro 



■ (B.5) 



m?=ro 



All diagrams appearing in Eqs. (B.4) and (B.5) can be calculated (at d = 3 
and q — 0) using Appendix A with mf = Tq. 

Following Ref. [6] we invert ^+(ro, Uq), Eq. (B.l), and get the function ro(^+, Uq). 
Subtracting Stq = ro — Tq, Eq. (6), from ro(C+, Uq) and letting ci — 3 we obtain 
Eq. (33). Below Tc we get ro(^-, uo), Eq. (23), at d = 3 as described in Ref. [6]. 
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